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ABSTRACT 

We have observed variations in the optical emission hnes from the X-ray binary Circi- 
nus X-1. These variations may be attributed both to orbital variations and to long 
term secular changes in line strength. We have detected double-peaked Ha emission 
lines on two occasions, providing the first direct evidence for an accretion disk in the 
system. The separation of the peaks was different on the two occasions, suggesting 
that the disk might have a different size. The equivalent width of the emission lines 
dropped by more than a factor of three between 1999 and 2000; this continues the 
trend seen in earlier data, so that the Ha equivalent width has now declined by a fac- 
tor of twenty since 1976. The emission lines do not appear to show signature of orbital 
motion, except for the data taken near phase 0, which show a significant velocity shift. 

We have observed an absorption component to the He i lines on one occasion. We 
suggest that, unlike the P Cygni profiles seen in X-ray spectra, this absorption does 
not arise in the accelerating zone of a radiatively driven wind. Instead, the absorption 
arises in material previously ejected from the system. It was only seen on this one 
occasion because the strength of the emission line had dropped dramatically. 
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1 INTRODUCTION 

Cir X-1 is a highly unusual X-ray binary, whose nature has 
been a puzzle for many years. Since the early 1970s, Cir X-1 
has shown erratic X-ray properties: its light curve differed 



of Cir X-1 with the supernova remnant suggests that the 
system may be a young (< 10^ y old) runa way system from 
a supernova explosion (Stewart et al. 1993). 

The optical counterpart to Cir X-1 was i dentified as a 



dramatfcally each Lime it was ubseived. Never Llieless, Llieie 



highl y-reddened star with strong Ha emission ( Whelan et al 



is a pe riodic modulation of 16.6 d (Kaluzicnaki et al. 1976 ^ 



which is believed to be the orbital period of the binary. The 
compact star in Cir X-1 is a neutron star, inferred from the 
Type I X-ray bursts observed in a brief episode (Tennant, 
Fabian & Shafer 1986). Further support is provided by the 
Z-source behaviour seen recently in X-ray timing data from 
the Rossi XTE satellite (Shirey, Bradt & Levine 1999). 

The radio counterpart of Cir X-1 is located 25' from the 
centre of the supernova remnant G321.9— 0.3, and i s appar- 
ently conn ected to the remnant by a radio nebula (Haynet 
et al. 1986). It shows flares at t he same period as the X-ray 



modulation ( Haynes et al. 197^) . Ci r X-1 has two arcminute- 
scale radio jets (^tewart et al. 199^), and an arcsecond-scale 



asymmetric jet s uggests the presenc e of relativistic outflow 
from the source (Fender et al. 1998). The close association 
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1977). This object was later shown to consist of three stars 
within a radi us of l"5, the southernmost of wh ich is the true 
counterpart ( Moncti 1992 ; Duncan et al. 1993| ). In 1997, we 

obtained spectra of Cir X-1 near apastronQ, using the Anglo- 
Australian Telescope (Johnston, Fender & Wu 1999, here- 
after Paper I). We detected an asymmetric Ha line which 
can be decomposed into two components, a narrow one and 
a broad component which is blue-shifted with respect to the 
narrow component. Comparison with archival spectra from 
the AAT and from HST showed that an asymmetric emission 
line has been present for the past twenty years. No previous 
observations have seen any evidence of an accretion disk. 

In this paper, we report on new optical spectroscopic 
observations of Cir X-1 obtained during 1999 and 2000. 



t Here and throughout, we use the words "periastron" and "apas- 
tron" to refer to phase and phase 0.5 rcsDcctivclv. as c alculated 
from the timing of the radio flares (Stewart et al. 1991). 
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Table 1. Journal of observations of Cir X-1. The columns show the UT date of observation, the Julian date of the midpoint of the 
observation, the telescope and instrument used, the mean phase of the observation, the total exposure time, and the wavelength range 
and resolution of the spectra, as measured from the arc lines. The phase was calculated according the ephemeris of Stewart et al. (1991). 



UT Date JD Instrument Phase t^xp Wavelength Resolution 

(s) range (A) (A) 



1999 Jul 11 


2451371.045 


2.3m H 


h DBS 


0.880 


18900 


5600- 


-7500 


2.7 


1999 Aug 20 


2451410.910 


2.3m H 


h DBS 


0.292 


5400 


5400- 


-7330 


2.1 


1999 Aug 21 


2451411.870 


2.3m H 


h DBS 


0.350 


1800 


5400- 


-7330 


2.1 


1999 Aug 22 


2451412.901 


2.3m H 


h DBS 


0.411 


7200 


5400- 


-7330 


2.1 


2000 May 16 


2451681.045 


AAT - 


h RGO 


0.622 


28800 


5660- 


-7000 


1.3 


2000 May 22 


2451687.063 


2.3m H 


h DBS 


0.985 


28800 


5570- 


-7000 


2.4 


2000 Jul 12 


2451738.096 


2.3m H 


h DBS 


0.073 


12600 


5580- 


-7510 


4.8 


2000 Jul 13 


2451739.013 


2.3m H 


h DBS 


0.129 


19800 


5508- 


-7510 


4.5 



2 OBSERVATIONS AND DATA REDUCTION 



3 RESULTS 



Cir X-1 was observed on multiple occasions during 1999 and 
2000, using the Double Beam Spectrograph (DBS) on the 
ANU 2.3-m telescope, and the RGO Spectrograph on the 
3.9-m AAT. A complete log of observations is given in Ta- 
ble |l[ 

The 2000 May 16 observations were carried out with 
the RGO spectrograph in combination with the MITLL3 
chip and a grating with 1200 grooves mm~^, resulting in 
a dispersion of 1.3 A pixel" ^. All other observations were 
carried out using the DBS with the standard SITe 1752 x 
532 CCD with a 600 grooves mm^^ grating. The very red 
colour of the source means that no significant signal was 
detected in the blue arm of the DBS. Series of 1800 s object 
exposures were interspersed with CuAr arc-lamp exposures. 
A slit width of 1.5 or 2 arcsec was used, and the slit was 
always oriented north-south so that both Cir X-1 and star 
2 of Moneti (1992) were in the slit. The spatial scale on the 
detector was 0'.'91 for the DBS and 0'.'48 for the RGO, which 
means the spectrum of Cir X-1 was confused with that of 
star 2; keeping the same slit orientation means that at least 
this contribution was constant. Because our previous work 
(Paper I) had shown no significant features associated with 
star 2, we have not attempted to remove its contribution 
from the spectrum. 

The phases of the observations are shown in T able m 
as calculated from the ephemeris of Stewart et al. ( 1991^ , 
based on the timing of the radio flares. It can be seen that, 
while there is moderately good coverage of a range of or- 
bital phases, there is a year separating the earliest spectra 
from the last spectra, so there are potential secular varia- 
tions which may confuse any features arising from the orbital 
variation. 

The IRAF software suite was used to remove the bias 
and pixel-to-pixel gain variations from each frame. As we 
had multiple consecutive observations of the same object, 
cosmic ray events were removed using the technique de- 
scribed by Croke (1995), as implemented in FIGARO. The 
wavelength calibration was performed in iraf by fitting a 
low-order polynomial to the arc line wavelengths as a func- 
tion of pixel number; the rms scatter of the fits was typically 
better than 0.1 pixel. 



The spectra all show a strong Ha emission line; the shape 
however changed significantly between the various observa- 
tions. We can see both long-term secular variation, in the 
equivalent widths of the lines, and variations which seem to 
be correlated with the orbital phase of the observations, in 
the line-profile morphology. Unfortunately, with sparse or- 
bital coverage, it is difficult to disentangle the latter from 
the former. Only with observations within a single orbit, 
and repeated observations at the various orbital phases, will 
we be able to be certain about some of the phenomena we 
describe here. However, with that caveat in mind, we will 
describe the data at hand. 



3.1 Ha line shape 

As was seen in previous observations (Paper I), the Ha 
line shows multiple components. In particular, the spectrum 
taken on 1999 July 11 shows a clear double-peaked profile. 
Such profiles are well known in the optical spectra of soft 
X-ray transients and cataclysmic variables, where they are 
usually take n to indicate the presence of an accretion disk 
( ^mak 1981 ) . A natural interpretation for the double-peaked 
lines in our spectra would therefore be that there is an ac- 
cretion disk present, although other possible interpretations 
are discussed below (§ [l.3| ). This is the first time evidence 
of the accretion disk has been seen in the optical spectra of 
Cir X-1. 

We fit Gaussian profi les to the Ha line, using the 
specfit package in iraf (Kriss 1994). After normalising 
each spectrum by a low-order polynomial fit to the contin- 
uum, we fit either two or three gaussians to the Ha line. The 
spectrum taken on 1999 July 11 clearly showed a double- 
peaked profile, in addition to the broad blue-shifted compo- 
nent, so we fit three gaussians to this spectrum. The reduced 
for 1/ = 85 degrees of freedom changes from = 4.9 for 
a two-gaussian fit to ~ 1-1 for three gaussians. The spec- 
trum from 2000 May 16 also required an extra narrow com- 
ponent to achieve a good fit (xS = 1-52 for a two-gaussian 



fit, Xv ~ 1-23 for three gaussians). The other spectra were 
fit with two gaussians. The details of the fits are listed in 
Table |. 

The data and the fits are shown in Figure |l|, with the 
spectra arranged in order of phase. Despite the spectra hav- 
ing been taken over a span of more than a year, there appears 
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Table 2. Fits to the emission line profiles. Column 2 shows the mean phase of the observation, column 3 shows the total equivalent width 
of the Ho line Wa,*! measured by direct summation of pixels in the emission line, and c olur nns 4 and 5 show the equivalent width of the 
He I A6678 and A7065 lines, measured from Gaussian fits to the line profiles (see Sect. |3.3| ). The He I lines in the 2000 May 16 spectra 
show P Cygni profiles; the equivalent widths listed are those of the emission component. Columns 6—11 show the results of Gaussian fits 
to the Ha line. The spectra for 1999 July 11 and 2000 May 16 were fit with three components, the other spectra with two. The fits are 
shown in Figure 1. 



W\ t W\ W\ Ha narrow component (s) Ha broad component 



UT Date 


Phase 


Ha 


A6678 


A7065 


velocity 


FWHM 


Wx 

(A) 


velocity 


FWHM 








(A) 


(A) 


(A) 


(kms"-'-) 


(kms-i) 


(kms"-'-) 


(kms-l) 


(A) 


1999 Jul 11 


0.880 


81 


3.5 ±0.4 


4.7 ±0.3 


-104 ± 10 
270 ± 10 


270 ± 16 
400 ± 15 


21 ±2 
40 ±1.5 


-600 ± 40 


550 ± 70 


9.1 ±1.1 


1999 Aug 20 


0.292 


80 


4.6 ±0.9 


6.7±1.0 


261 ±8 


490 ± 20 


47 ±4 


-90 ±60 


1010 ±60 


27 ±4 


1999 Aug 21 


0.350 


79 


1.7±1.0 


4.4 ±0.8 


350 ± 16 


290 ± 40 


24 ±6 


180 ± 30 


675 ± 50 


51±6 


1999 Aug 22 


0.411 


83 


3.2 ±0.8 


3.9 ±0.6 


330 ± 10 


320 ± 24 


30 ±5 


1040 ± 30 


690 ± 30 


49 ±5 


2000 May 16 


0.622 


12.2 


0.28 ±0.08 


0.66 ±0.09 


-30 ±5 
200 ±5 


210 ±40 
170 ± 30 


2.7 ±0.7 
3.1±0.7 


-420 ± 30 


630 ± 45 


6.6 ±0.9 


2000 May 22 


0.985 


10.1 


1.4±0.2 


0.3 ±0.1 


87 ±15 


270 ± 60 


1.5±0.5 


44 ±18 


1000 ± 50 


9.5 ±0.5 


2000 Jul 12 


0.073 


25.0 


2.0±0.5 


1.9 ±0.3 


307 ± 80 


580 ± 60 


11 ±15 


-20 ±400 


740 ± 230 


10 ±15 


2000 Jul 13 


0.129 


26.3 


2.6 ±0.5 


2.3 ±0.4 


290 ± 40 


440 ± 40 


13.8±3 


-115 ±70 


460 ± 80 


8.0 ±3 
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Figure 1. Line profiles of Ha, showing the Gaussian fits to the 
lines and their sum. The spectra have been normalised by a poly- 
nomial fit to the continuum. The spectra are shown in order of 
phase, with the spectrum taken at phase (1999 May 22) re- 
peated at phase 1 for clarity. The spectrum at phase is sym- 
metric, with a broad component on the blue wing appearing at 
phases 0.1—0.5. At phase 0.6 the line appears to be double-peaked 
(or flat-topped?), while clear double peaks are seen at phase 0.9. 



to be a clear trend of line shape with orbital phase. The spec- 
trum taken at phase (2000 May 22) is the only completely 
symmetric profile. At small phases, the line profile begins to 
show asymmetry in the line. Sometime after phase 0.4, the 
line shows a second peak in addition to the broad compo- 
nent (2000 May 16), which persists for the second half of the 
orbit (1999 July 11). We have fit two narrow components to 
the line at phase 0, although it is possible that the profile 
is fiat-topped instead of intrinsically doubl e-pe aked. We will 



discuss the implications of this in Section 4.3 



In Paper I we presented five archival spectra from the 
AAT and one from HST; we can re-examine these spectra 
for evidence of the same pattern in the line profiles. Three 
spectra (1976 May, 1978 August and 1997 June) were taken 
near apastron, at phases 0.54, 0.52 and 0.51; none of them 
shows a double-peaked profile. The spectra taken in the first 
half of the orbit all show an asymmetry on the blue wing; 
no spectra were taken during the second half of the orbit 
(phases 0.6-0.9), so we do not know whether the double- 
peaked line has been consistently present. The spectrum 
taken near periastron (1977 September, phase 0.95) does not 
show a symmetric profile; the line is still asymmetric on the 
blue wing. Thus we can say very little about the persistence 
of the pattern we find in our spectra from 1999-2000. 



3.2 Secular variation in Ha equivalent width 

The equivalent width Wx of the Ha line is shown in Ta- 
ble ^ and its evolution is plotted in Figure ^ The equiva- 
lent width dropped dramatically between 1999 (MJD 51370- 
51412) and 2000 (MJD 51680-51738), from 80 A to 25 A, 
and the drop is independent of the phase of the observa- 
tions. This drop coincides with a drop in the X-ray fiux 
observed by RXTE, where the baseline fiux between fiares 
(~ 100 cts-i for the All-Sky Monitor (ASM) on RXTE) be- 
gan a slow decline around 2000 March (MJD 51600). This 
downward trend of the Ha equivalent width is consistent 
with the trend seen over the last 25 years (Paper I, Table 2), 
where the observations in 1976 May found W\ = 580 A, de- 
clining to 120 A in 1997. Rough fiux-calibration of this data 
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Figure 2. Variation of tlie X-ray flux and equivalent width of the 
Ha line with time. The top section of each panel shows the 1-day 
average ASM flux measured by XTE, the lower section shows the 
equivalent width of the Ha line. The entire RXTE/ASM light 
curve for Cir X-1 is shown in the top panel, while an expanded 
version of the light curve since 1999 is shown in the bottom panel. 
The "baseline" X-ray flux, between flares and dips, began to de- 
cline from a value of about 1.2 Crab (shown by the dashed line) to 
0.7 Crab by 2000 July (MJD 51740). The Ha equivalent widths 
for 1995 (phase 0.18) and 1997 (phase 0.51) are from Paper I. 
The drop is independent of phase. 



and the archival data shows that this change in the equiva- 
lent width is not due to an increase in the brightness of the 
continuum, and therefore must reflect a real and dramatic 
decrease in the brightness of the emission line. 



3.3 He I lines 

The spectra all show emission lines of He I A6678 and A7065. 
In most cases, these lines are quite faint; the equivalent 
width of the He I lines also dropped dramatically between 
1999 and 2000. We measured the equivalent width of the 
lines by fitting Gaussians in the same manner as for the Ha 
lines. For most of the spectra, a single gaussian for each line 
was sufficient. For the 1999 July 11 lines, the He I lines were 
double peaked just like the Ha line. The 1999 August 20 
line required two components for a satisfactory fit. 

The spectrum from 2000 May 16 clearly showed P Cygni 
profiles in both the He I lines (Figure P Cygni profiles 
were recently discovered in X-ray spe ctra of Cir X-1 taken 
by Chandra ( Brandt fc Schulz 2000| ). The peak-to- valley 
distance of the profile in the optical spectra is 7.62 A, or 
340 kms"^ This is a much smaller velocity than that seen 
in the X-ray spectra (2000-3000 kms"^). 

Figure ^ shows the equivalent width of the He I lines 
as a function of phase. In order to remove the dramatic 



99 Aug 20 



00 May 1 6 




6660 6580 6700 
Wavelength (A) 



7060 7080 
Wavelength (A) 



Figure 3. Selected line proflles of He I, showing the flts to the 
lines. The spectra have been normalised by a polynomial fit to 
the continuum; the wavelengths of the flts were constrained to 
have the same velocity and width. The left side of the flgure 
shows the A6678 line, the right side the A7065 line. The 1999 
July 11 spectrum shows a double-peaked profile, the 1999 August 
20 spectrum shows an asymmetric line, and the 2000 May 16 
spectrum shows a P Cygni profile. 



secular variation which occurred between 1999 and 2000, 
we normalised the He I equivalent widths by the equivalent 
width of the Ha line (Table ^, column 3) . The results show 
a smooth variation with phase; both transitions are of sim- 
ilar strength, except for the abrupt reversal in the phase 
spectrum, when the A7065 line becomes significantly weaker 
than the A6678 line. 



3.4 Velocities of the emission lines 

The measured velocities of the emission line do not seem 
to conform to any sensible pattern. In particular, there is 
no discernible signature of orbital motion over most of the 
orbit. This implies that the site of the emission is changing 
throughout the orbit, and possibly over the one year span- 
ning our observations - see Section ^ for a discussion of the 
implications for a model for the system. 

However, we do see significant orbital motion in the 
data taken near periastron. We divided the data taken on 
2000 May 22 into four equal segments of 7200 s duration and 
performed the same analysis on each segment. The results 
are shown in Fig. |^ Both components of the Ha line, as 
well as the He I A6678 line (not plotted), show an increase in 
velocity of 120 kms~^ over six hours. This certainly implies 
that at this point in the orbit the emission line is tracking 
the orbital motion; in a highly elliptical orbit, nearly all the 



veloc ity swing happens right near periastron ( Tauris et al 



199g) 



With only four velocity points, we cannot place any 
interesting constraints on the orbital parameters. Two rep- 
resentative curves are plotted with the observed points; it 
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Figure 4. Equivalent width of the He I Hnes, normalised by the 
Hq equivalent width, as a function of phase. The filled squares 
show He I A6678, the crosses Hel A7065. The dates of the obser- 
vations are indicated above the plot. The spectrum at phase 0.62 
(2000 May 16) showed P Cygni profiles in the Hel lines: only 
the equivalent width of the emission component has been plotted 
here, thus underestimating the intrinsic line strength. 
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may22c 
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can be seen that there is essentially no constraint on the 
orbital inclination, though the eccentricity must be high to 
produce the observed sharp swing in velocity. 



3.5 Timing of periastron 

It is important to note that the phases have been calculated 
using the ephemeris of Stewart et al. (1991). Phase refers 



to the onset of radio flares, and how this relates to the X-ray 
dips observed by RXTE is not clear. We have examined the 
quick- look results provided by the RXTE/ ASM team to see 
if we can compare the observed times of the X-ray dips with 
the predicted times from the Stewart et al. ephemeris. The 
periastron passage on 2000 May 22 did not have good ASM 
coverage; both the preceding and the following X-ray dips 
appear to have occurred 0.4-0.5 d later than the phase 
predicted by the Stewart et al. ephemeris, corresponding to 
a shift of 0.024-0.03 in phase. Given that the ephemeris was 
based on radio flares, and that the X-ray light curve changes 
so dramatically, this agreement is reasonably good. In future 
work we hope to further investigate the exact relationship 
between the X-ray, radio and optical behaviour. 



Figure 5. Velocities of the fitted components near phase 0. The 
data from 2000 May 22 has been divided into four spectra of 
duration 7200 s; the rapid swing of velocity in both the narrow 
(solid line) and broad (dashed line) components is evident in the 
expanded section to the lower right. The dotted line shows the 
second narrow (blueward) peak in the double-peaked spectra on 
99 July 11 and 2000 May 16. The dotted and dot-dashed lines 
show two representative orbital curves, assuming Mi = 1.5 Mq, 
M2 = 3.0 Mq; the dotted line has i = 45deg, e = 0.9, while the 
dot-dashed line has i = 85deg, e = 0.87; both curves have the 
longitude of periastron set to 210 deg and a systemic velocity of 
— 120 kms~^ and —220 kms"'^ respectively. 



strong X-ray variability after phase 0. Between phases 0.5- 
0.9, the accretion disk becomes Keplerian and there is steady 
accretion. The disk is then disrupted by tidal forces during 
the following periastron passage. 

Based on the consistently high velocities seen in all the 
data up to that time, we suggested in Paper I that Cir X-1 
might have a very high radial velocity. Our new spectro- 
scopic data do not show evidence of high line-of-sight ve- 
locities, though the argument about the velocity from the 
association with G321.9— 0.3 is still valid. 



4 DISCUSSION 

4.1 A model for the system 

In Paper I, we proposed a theoretical model for Cir X-1. 
The binary is a low-mass neutron-star binary with an ultra- 
eccentric binary orbit (e ^ 0.7, possibly as high as 0.9; 
Tauris et al. 1999). During periastron passage, the com- 
panion star (which is a subgiant of about 3-5 Mq) over- 
fills its Roche-lobe, causing a transfer of mass at a super- 
Eddington rate, which in turn drives a strong matter out- 
flow. After periastron, mass transfer from the companion 
ceases, but accretion continues at a near-Eddington rate as 
the neutron star captures the residual matter in its Roche- 
lobe. An accretion disk gradually forms, which is initially 
non- Keplerian and geometrically thick. This change between 
quasi-spherical and disk accretion causes the change between 



4.2 Sites of emission 

Our data and the archival data show a persistent broad com- 
ponent in the Ha line, indicating the presence of a high- 
velocity flow. One possibility for the line-formation region is 
the accretion flow close to the neutron star. The X-ray lumi- 
nosity of Cir X-1 reaches the Eddington limit at the phases 
shortly after periastron. This implies a mass accretion rate 
M Z2x 10^** gs~\ The quantity n^R^ is given by 

n,R^ « l.OxlO^M ^„ A 

\2 X 10'**gs-V 

A000 kms"^\ _i _i 
X ( I X cm , (1) 

where Uc is the electron number density, x parametrises the 
degree of spherical accretion, and its value is unity for spher- 
ical accretion. Hence the ionisation parameter 
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(in cgs units) 



(2) 



implying a temperature T ~ 10^ K for the accreting matter 
(see Hatchett et al. 1976). The strong X-ray emission from 
the source would therefore completely ionise the accretion 
gas with inflowing velocity with v ~ 1000 kms^^, making it 
unlikely to emit H Balmer lines. 

Alternatively, the broad component of the Ha line could 
be emitted from the outflow matter. A characteristic of a ra- 
diatively driven outflow is that the outflowing matter is ac- 
celerated from a slow initial velocity and eventually reaches 
a terminal velocity. As there is a negative gradient from the 
observer to the outflow matter, the blue wing of the line 
suffers very strong self-absorption. The line would therefore 
show a P Cygni profile. Since our optical data do not show 
the P Cygni profile typical of such a radiatively driven out- 
flow (with the exception of one set of spectra: see below), the 
gradient of the line-of-sight velocity in the line formation re- 
gion must be positive, i.e. the outflow is decelerating. If the 
broad component originates from a high-velocity outflow, 
the line formation region is in the decelerating zone, and 
the accelerating zone is behind the photosphere. Such out- 
flow might be caused by an ejection of high- velocity material 
during the violent mass transfer epoch at the periastron. 

On one occasion, we do see an absorption component in 
the blue wing of the He I lines, though not in the Ha line. We 
suggest that, instead of arising in a radiatively driven out- 
flow, we are seeing absorption from the outflowing material 
far from the system, superimposed on emission lines from 
the accretion disk. The absorbing matter consists of mate- 
rial ejected from the system in previous episodes of violent 
mass transfer, and the absorption is occurring at some dis- 
tance from the system, where the outflow has decelerated 



signific antly from th e 2000 3000 kms se e n in th e X-ray 

spectra! 



The absorption is seen only in spectra from 2000 May 
16, which is when the He I emission from the disk was weak- 
est (Table ^). We suggest that the absorption in the outflow 
is masked in the other spectra by the strength of the emis- 
sion line, which is also why we do not see similar absorption 
in the Hq line. 



4.3 E irot dotoction of an accretion diak? 



Our detection of the double-peaked lines formed in an accre- 
tion disk is strong evidence for the eccentric binary model for 
Cir X-1. The double-peaked lines were seen twice, in spectra 
taken almost a year apart, at phases 0.88 and 0.62. Spectra 
taken in the first half of the orbit all show an asymmetric 
Ha emission line, while the spectrum taken at periastron 
(phase 0.98) shows a completely symmetric line. 

Double-peaked lines are often found in cataclysmic vari- 
ables (e.g. IP Peg, see Marsh 1988), which are binaries with a 
white dwarf accreting material from a low-mass companion 
star, and in black-hole binaries during the high-soft state 
(e.g. GX339-4, see Soria, Wu & Johnston 1999) or qui- 
escence (e.g. A0620-00, Johnston, Kulkarni & Oke 1989). 
They are also found in the accretion-disk corona sources 



rare amongst other X-ray binaries, and X-ray pulsars do not 
typically show double-peaked optical line^. 

It was shown (Wu et al. 2001) that double-peaked lines 



can be formed in a temperature-inversion layer on an opaque 
accretion disk when the disk is irradiated by soft X-rays. The 
X-ray spectrum of Cir X-1 is soft, dominated by a black- 
body component w ith an effective temperature T ~ 2 keV 
( ^hirey et al. 199'9b ). Thus the double-peaked lines observed 
in Cir X-1 probably indicate the presence of an accretion 
disk. Moreover, a significant part of the disk surface is not 
obscured from the central X-ray source so that the irradia- 
tively induced temperature inversion can occur. 

It is interesting to note that the 1999 July (Fig. |l|, 
phase 0.88) observation, which was the first spectrum to 
show double-peaked lines, was taken in the cycle after an 
orbit during which Cir X-1 did not have a post-periastron 
X-ray flare; only a handful of orbits observed by XTE had 
not shown such a flare before then. It is possible that this 
has some bearing on the visibility of the optical lines from 
the accretion disk, for instance if the absence of an X-ray 
outburst means that matter which usually obscures the disk 
was more transparent than usual. Alternatively, the accre- 
tion rates during this observational epoch were low enough 
to allow a mature Keplerian disk to develop within an orbital 
cycle. 

The double-peaked lines seen in 2000 May (Fig. |l|, phase 
0.62) did not occur after such an anomalous orbit. However, 
by that stage the X- ray flux observed by RXTE had begun 
to drop (see Sec |3!|, Fig. |), which might also imply that 
veiling material was thinner. 

It is possible that, instead of being an intrinsically 
double-peaked line, the profile observed on 2000 May 16 
is actually flat-topped. Such lines have seen in several black 



hole candid ates, such as GX 339—4 ( Smith et al. 199 i: ; 3oria 
et al. 1999|) and GRO J 1655 -40 ( |Soria et al. 2000| ), where 



they are inter preted as arising in a hor izontal wind launched 
from a disk (Murray & Chiang 1997). Possibly the line in 
Cir X-1 arises in a similar fashion; however, this requires 
the emission site to change from a spherical outflow (phases 
< 0.5; Fig. |l|) to a horizontal wind (phase 0.6) to a disk 
(phase 0.88). The interpretation of the line as a barely- 
resolved double peak has at least the virtue of being a sim- 
pler explanation. 

Assuming the double-peaked lines arise in an accretion 
disk, then the separation reflects the velocity at the outer 
edge of the emission region of the accretion disk (3mak 



1981). The separation of the peaks was very different on 



the two occasions the double-peaked profile was observed 
(370 kms"^ in 1999 July, 230 km s^^ in 2000 May), which 
suggests that the accretion disk had a very different size. 
This could be due to the different phase of the observation, 
implying the disk is shrinking as periastron approaches. Fig- 
ure o shows a comparison between the observed peak sepa- 
rations and a model where the size of the disk is truncated 
by the size of the Roche lobe in an eccentric orbit. This rep- 
resents the largest possible size for the accretion disk, and 
hence a lower limit to the separation of the two peaks. The 
fact that the observed velocities are above this limit implies 



(e.g. 2 4 1822-371, Mason et al. 1982). They are, however. 



^ Though note th at 2A 1822—371 was recently shown (Jonker . 
van der Klis 2001 ) to be an X-ray pulsar! 
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Figure 6. Sketch of the separation of the double-peaks as a func- 
tion of orbital phase, for some representative orbital parameters, 
compared with our observed values. The top panel shows the 
separation between the two stars (top line) and the outer radius 
of the disk, where we assume the disk is limited by the size of 
the Roche lobe of the compact star and its outer radius is trun- 
cated at 0.7 Rij. The masses of the two stars are Mi = 1.5 Mq, 
M2 = 3.0 Mq, and the orbital eccentricity and inclination are 
e = 0.9 (solid lines) and e = 0.7 (dotted lines). The bottom panel 
shows the predicted peak separation for such a disk, for a range 
of orbital inclinations: (from the bottom) i = 30 deg, 45 deg and 
75deg. The open squares show the measured peak separation for 
the two dates on which they were observed. 



that either the disk has not reached the truncation radius, 
or that the emission region for Ha is not at the outer edge 
of the disk. The former would imply that, after the disk 
is disrupted at periastron, it has not had time to grow to 
the truncation radius before the Roche lobe begins to con- 
tract again (Fig. P), which would imply that the viscous 
timescale for the growth of the disk is longer than the dy- 
namical timescale, which is 8 d for half an orbit. 

Alternatively, the change in peak separation could be 
the result of a secular variation in the size of the emitting 
region; or the temperature of the disk could be varying, so 
that the location of the Ha emission site is changing. In this 
case, the smaller peak separation implies a hotter disk, with 
the emission site having moved towards the outside of the 
disk. 

Other explanations for the double-peaked profiles are 
possible which do not require the presence of an accretion 
disk. They could, for instance, arise in a bipolar jet, since 
a relativistic jet has bee n inferred from high -resolution ra- 
dio maps of the source (Fender et al. 1998). However, in 



SS 433, which has a jet with velocity ~ 0.3c, the separatio n 
of the components is hundreds of angstroms (Margon 1984), 
compared to the 5-8 A seen in our spectra. Such a small 
separation could only be produced by a relativistic jet if 
the jet was almost in the plane of the sky (in which case the 
transverse Doppler shift would cause both components to be 
redshifted). Thus any jet would need to be non-relativistic 
to produce the observed lines. 

Variable double-peake d lines have also be en seen in in- 
frared spectra of Cyg X-3 (Fender et al. 1999), where they 
were interpreted as arising in a disk-like wind outside the 
binary orbit. Such a disk-wind requires a large angular mo- 



mentum, which would be unlikely in a system like Cir X- 
1 with such a large orbital period. The double-peaks could 
arise from emission in different parts of the system; however, 
we note that the standard sites for such localised emission, 
such as the hot spot where the accretion stream impacts the 
disk, still require the presence of a disk in the system. 

Finally, it is possible that the central dip in the double- 
peaked profile is formed by absorption in front between us 
and the emitter, especially given the fact that we s ee a bsorp- 
tion in the He I lines on a different occasion (see § 3.3). Such 



lines are seen in symbiotic stars, such as CH Cygni, which 
shows t wo symmetric peaks s eparated by a deep central re- 
versal (Anderson et al. 1980); in this object, the profile is 
assumed to arise in a self-absorbed wind. However, in the 
case of the profiles seen in Cir X-1, the conditions for such 
absorption would be hard to produce. Since the profile is 
symmetric, the absorbers would need to be at a very similar 
velocity to the emitters; and since the central dip is broad 
and does not reach the continuum, the absorbers would need 
a wide velocity dispersion and only a small optical depth. 
The most probable interpretation of the profile is that it 
arises from the Doppler motions in an accretion disk. 



4.4 Long-term changes in the emission lines 

The brightness of the emission lines has decreased dramat- 
ically over the past 25 years, from W\ = 580 A in 1976 to 
25 A in 2000. Our interpretation is that the contribution of 
different emission sites to the emission lines has been shift- 
ing over the years. The contribution from the accretion disk 
is likely to have been essentially constant, since the temper- 
ature of the disk is limited by the accretion rate (assumed 
to be at Eddington) and the size of the disk is limited by 
the size of the orbit. The contribution to the emission line 
from other regions - the outflow, the companion star - can 
change dramatically depending on the state of the system. It 
is possible that we have seen the disk for the first time pre- 
cisely because the contributions from other sites in the orbit 
have lessened, so the Ha line is now dominated by emission 
from the accretion disk. 
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